Dilepton production in proton-proton and Pb+Pb collisions at snn = 2.76 TeV 
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We study e + e~ pair production in proton-proton and central Pb+Pb collisions at y/swN = 
2.76 TeV within two models: an extended statistical hadronization model (SHM) and the Parton- 
Hadron-String Dynamics (PHSD) transport approach. We find that the PHSD calculations roughly 
agree with the dilepton spectrum from hadronic sources with the 'cocktail' estimates from the sta- 
tistical hadronization model matched to available data at LHC energies. The dynamical simulations 
within the PHSD show a moderate increase of the low mass dilepton yield essentially due to the 
in-medium modification of the p-meson. Furthermore, pronounced traces of the partonic degrees of 
freedom are found in the PHSD results in the intermediate mass regime. The dilepton production 
from the strongly interacting quark-gluon plasma (sQGP) exceeds that from the semi-leptonic de- 
cays of open charm and bottom mesons. Additionally, we observe that a transverse momentum cut 
of 1 GeV/c further suppresses the relative contribution of the heavy meson decays to the dilepton 
yield, such that the sQGP radiation strongly dominates the spectrum for masses from 1 to 3 GeV, 
allowing a closer look at the electromagnetic emissivity of the partonic plasma in the early phase of 
Pb+Pb collisions. 

PACS numbers: 25.75.-q, 13.60.Le, 14.40.Lb, 14.65.Dw 



I. INTRODUCTION 

Dilcptons, i.e. correlated electron and positron or 
pairs, are one of the key observables in experiments 
for ultra-rclativistic nuclear collisions since they are emit- 
ted during the whole evolution of a collision and interact 
only electromagnetically and thus very weakly with the 
strongly interacting partonic or hadronic medium created 
in the collisions. Also, dileptons of different invariant 
masses are dominantly produced from different stages of 
a relativistic nuclear collision, and this provides the pos- 
sibility to probe the properties of the produced hot and 
dense matter at various conditions by measuring the dif- 
ferential dilepton spectra. 

The invariant mass spectrum of dilcptons can be 
roughly divided into 3 different regions which are dom- 
inated by different physics. In the low mass region 
(M e + e - < 1 GeV), the radiation is dominated by 
the decays of light mesons (consisting of u, d and s 
(anti)quarks), while in the intermediate mass region (1 
GeV < M e + e - < 3 GeV), the dominant hadronic con- 
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tribution to the invariant mass spectrum stems from the 
decays of open charm mesons. Above the J /tp peak, the 
dilepton spectrum is first dominated by open beauty de- 
cays and then by the initial state Drcll-Yan radiation. 
Besides these sources of dileptons, the radiation from the 
strongly interacting Quark-Gluon-Plasma (sQGP) [1] as 
well as some other more exotic sources like simultaneous 
interactions of four pions [2-5] can also give significant 
contributions, particularly in the intermediate mass re- 
gion. These partonic and hadronic channels have been 
studied in detail in Refs. [6, 7] at top Super-Proton- 
Synchrotron (SPS) and Relativistic-Heavy-Ion-Collidcr 
(RHIC) energies, and it has been found that the par- 
tonic channels clearly dominate over multi-pion sources 
in the intermediate dilepton mass regime. 

In the present study, we include all known (lead- 
ing) dilepton sources to study and compare the mag- 
nitude of the radiation from the QGP with that from 
other (conventional) sources in heavy ion collisions at the 
Large-Hadron-Collider (LHC) energies. Since at present 
there is no single model that could address reliably all 
of the above-mentioned sources for dilpeton production, 
we thus employ in this analysis different approaches to 
evaluate the invariant mass spectrum from the different 
sources. Specifically, we concentrate on two models (and 
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their extensions): The Parton-Hadron-String-Dynamics 
(PHSD) approach [8, 9] and the extended Statistical 
Hadronization Model (SHM) [10]. The PHSD model is a 
relativistic transport approach developed, tested and well 
suited for studying dynamical partonic and hadronic sys- 
tems in heavy ion collisions from low to ultra-relativistic 
energies of ^/snn=200 GcV, while the SHM is a statis- 
tical hadronic model suitable for describing in detail the 
relative yields of final state hadrons. 

In the case of proton-proton (pp) collisions, they are 
dominated by the creation and decays of various hadrons, 
and the SHM should be a suitable model for estimating 
the hadronic freeze-out 'cocktail' contribution to dilep- 
tons. Indeed, this was confirmed by our simulations 
in Ref. [10] for p+p collisions at ^/snn= 200 GeV. On 
the contrary, in heavy-ion collisions a rather long-living 
strongly interacting dynamical fireball is formed and 
dilcptons arc emitted from the created charges over an ex- 
tended period of time. For such systems the SHM might 
not provide a good description because the properties, i.e. 
the masses and widths, of hadrons might change in the 
hot and dense nuclear environment as a function of time, 
which can not be taken into account in the SHM. Ad- 
ditionally, the dilepton signal from the QGP during the 
early reaction phase is expected to outshine the back- 
ground radiation due to the cocktail of hadron decays. 
For this purpose, we will employ the PHSD transport 
approach which can account for the medium-dependent 
properties of partons and hadrons and their multiple in- 
teractions as well as for the dynamical evolution of the 
system in general, including the transition to the partonic 
phase and the radiation from the sQGP. 

Wc note that a direct comparison between results from 
the PHSD and the SHM has shown that the low mass 
(M < 1 GeV) spectra of e + e~ -pairs in Au+Au collisions 
at v^nn~ 200 GeV deviate from each other by only up 
to 20% [7, 10], i.e. either the time-evolution of the system 
plays only a minor role in the dilepton emission or the 
subsequent hadronic evolution to a large extent washes 
out the details of the dynamics in the dilepton invariant 
mass spectrum. We will study this issue again in detail in 
the first part of present study and compare explicitly our 
results from both the static and dynamical model in or- 
der to constrain the dilepton yield from known hadronic 
sources. We will then proceed to our main goal, i.e. the 
identification of the QGP signal from the spectrum of 
dilcptons produced in Pb + Pb collisions in the LHC en- 
ergy regime. 

In the following Section (Sect. II), we will briefly re- 
call the main concepts of the PHSD transport approach 
and list the partonic and hadronic sources of dileptons in 
PHSD. Next, we describe the SHM and the way hadron 
yields and spectra are evaluated in Sect. III. Our treat- 
ment of D and B meson production and energy loss is 
described in Sect. IV, whereas the modeling of charmo- 
nia production in elementary and in nucleus-nucleus col- 
lisions is outlined in Sect.V. Results from our study are 
presented in Sect. VI. We first include a discussion of 



dilepton production in the baseline p+p collisions. Wc 
then continue with results from the SHM as well as from 
the dynamical calculations of Pb+Pb collisions within 
the PHSD. Furthermore, the effect of in-medium modifi- 
cations of vector mesons on the low-mass dilepton spec- 
trum at y / SNN=2.76 TcV is investigated. Finally, the 
influence of a transverse momentum cut on the dilepton 
spectrum as a possible way to enhance the QGP signal 
is explored. We conclude our study in Sect. VII with a 
summary of our findings and a discussion of open prob- 
lems. 



II. THE PARTON-HADRON-STRING 
DYNAMICS (PHSD) TRANSPORT APPROACH 

The PHSD model [8, 9] is an off-shell transport model 
that consistently describes the full evolution of a relativis- 
tic heavy-ion collision from the initial hard scatterings 
and string formation through the dynamical deconfine- 
ment phase transition to the quark-gluon plasma as well 
as hadronization and to the subsequent interactions in 
the hadronic phase. In the hadronic sector, the PHSD is 
equivalent to the Hadron-String-Dynamics (HSD) trans- 
port approach [11-13] that has been used for the descrip- 
tion of pA and AA collisions from SIS to RHIC energies 
and has led to a fair reproduction of measured hadron 
abundances, rapidity distributions and transverse mo- 
mentum spectra. In particular, as in the HSD, the PHSD 
incorporates off-shell dynamics for vector mesons [14] and 
a set of vector-meson spectral functions [15] that covers 
possible scenarios for their in-medium modifications. In 
PHSD the transition from the partonic to hadronic de- 
grees of freedom is described by covariant transition rates 
for the fusion of quark-antiquark pairs to mesonic reso- 
nances or three quarks (antiquarks) to baryonic states, 
i.e., by the dynamical hadronization [16]. Note that due 
to the off-shell nature of both partons and hadrons, the 
hadronization process obeys all conservation laws (i.e., 
the 4-momentum conservation and the flavor current con- 
servation) in each event, the detailed balance relations, 
and the increase in the total entropy S. The transport 
theoretical description of quarks and gluons in the PHSD 
is based on a Dynamical QuasiParticlc Model (DQPM) 
for partons that is constructed to reproduce the lattice 
QCD (1QCD) results for a quark-gluon plasma in ther- 
modynamic equilibrium. The DQPM provides the mean- 
fields for gluons/quarks and their effective 2-body inter- 
actions that are implemented in the PHSD. For details 
about the DQPM model and the off-shell transport, wc 
refer the reader to the review in Ref. [17]. 

We stress that a non-vanishing width in the partonic 
spectral functions is the main difference between the 
DQPM and conventional quasiparticle models [18]. Its 
influence on the collision dynamics can be seen in the 
correlation functions, which, in the stationary limit in- 
volve the off-diagonal elements of the energy-momentum 
tensor T kl and thus give rise to the shear viscosity r\ of 
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FIG. 1: (Color on-line) Pseudo-rapidity distribution of 
charged hadrons at midrapidity as a function of the number 
of participants N par t from PHSD (solid line) in comparison 
to the data from the ALICE Collaboration [23] for Pb+Pb at 
2.76 TeV. 



the medium [19]. A sizable width is then essential for 
obtaining a small ratio of the shear viscosity to entropy 
density rj/s, which results in a roughly hydrodynamical 
evolution of the partonic system in PHSD [16]. The two- 
particle correlations resulting from the finite width of the 
parton spectral functions are taken into account dynam- 
ically in the PHSD by means of the generalized off-shell 
transport equations [14] that go beyond the mean field 
or Boltzmann approximation [17, 20]. 

We recall that the PHSD approach has been tested 
from low SPS to top RHIC energies against the measured 
rapidity spectra of various particle species, transverse 
mass distributions, differential elliptic flow of charged 
hadrons as well as their quark-number scaling [21]. More 
recently, higher harmonics in the azimuthal distribution 
of charged hadrons in the plane perpendicular to the 
beam direction has also been examined [22]. The de- 
scription of the various data sets has been found to be 
surprisingly good for all bulk observables. 

For nucleus-nucleus collisions at LHC energies, one 
might wonder if some new features, which are not prop- 
erly described by the PHSD approach within the stan- 
dard settings of Ref. [9], could appear since in the ini- 
tial state a color glass condensate might become impor- 
tant [24] or the nuclear shadowing could be different 
from the extrapolations given in Ref. [9], etc. To shed 
some light on these questions, we have tested the default 
PHSD version by using it to study Pb+Pb collisions at 
^/snn = 2.76 TeV and compared the results to data from 
the ALICE Collaboration [23]. In particular, we have ex- 
amined the pseudo-rapidity distribution of charged par- 
ticles as a function of centrality. In Fig. 1 we display 
the results for dN c /dr) at midrapidity from the default 
PHSD calculations in comparison to the ALICE data as 
a function of the number of participants N part that has 



been determined dynamically in the PHSD calculations. 
Quite acceptable agreement is seen, indicating that the 
bulk parton dynamics is not much different at top RHIC 
and LHC energies. A similar observation has also been 
made in Ref. [25] where the UrQMD transport approach 
(without explicit partonic degrees of freedom) was shown 
to give a reasonable reproduction of this observable as a 
function of centrality. Thus we continue with electromag- 
netic probes using the default PHSD version of Ref. [9] 
without readjusting any PYTHIA parameters used in the 
PHSD for the initial state hard scattering processes. 



A. Partonic sources of dileptons in PHSD 

The PHSD approach so far has been employed for 
dilcpton production from pp to Au + Au collisions at 
SPS [6] and RHIC energies [7]. As in Refs. [6, 7] dilep- 
ton radiation from the constituents of the strongly in- 
teracting QGP proceeds via the following elementary 
processes: the basic Born q + q annihilation mecha- 
nism, gluon Compton scattering (q + g — > 7* + q and 
q + g — > 7* + q), and quark and anti-quark annihila- 
tion with the gluon Bremsstrahlung in the final state 
(q + q — > g + 7*). In the on-shcll approximation, one 
would use pcrturbative QCD cross sections for the pro- 
cesses listed above. However, in the strongly interacting 
QGP the gluon and quark propagators differ significantly 
from the non-interacting propagators. Accordingly, we 
have calculated in Refs. [26, 27] the off-shell cross sections 
for dilcpton production in the partonic channels by off- 
shell partons, using the phenomenological parametriza- 
tions from the DQPM for the quark and gluon propaga- 
tors and their interaction strength. 

We have implemented the cross sections obtained in 
Refs. [26, 27] into the PHSD transport approach in the 
following way [7]: Whenever quark-antiquark, quark- 
gluon and antiquark-gluon collisions occur in the course 
of the Monte-Carlo simulation of the partonic phase in 
the PHSD, a dilcpton pair can be produced according to 
the off-shell cross sections [27]. In addition to the virtu- 
alities of the partons involved, the latter also depends on 
the energy density in the local cell, where the collision 
takes place, as according to the DQPM it governs the 
widths of the quark and gluon spectral functions as well 
as the strong coupling. 



B. Hadronic sources of dileptons in PHSD 

In the hadronic sector, the PHSD model is equivalent 
to the HSD transport approach [11-13]. The implemen- 
tation of the hadronic decays into dileptons (tt-, 77-, rf-, 
uj-, A-, eti-Dalitz, p — > , uj — > and </> — > l + l~) in 

HSD (and PHSD) is described in detail in Refs. [15, 28]. 
In contrast to the HSD approach - without explicit par- 
tonic degrees of freedom - the mesons in the PHSD are 
produced through the dynamical hadronization from the 
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partonic state. On the other hand, the subsequent inter- 
actions of the produced hadrons as well as the evolution 
of the hadronic 'corona' proceed in the same way as in 
the HSD. 

The PHSD off-shell transport approach is particularly 
suitable for investigating the different scenarios for the 
modification of vector mesons in a hot and dense medium. 
As in the HSD model, the PHSD approach incorporates 
the off-shell propagation of vector mesons as described 
in Ref. [14]. In the off-shell transport, the hadron spec- 
tral functions change dynamically during the propaga- 
tion through the medium and evolve towards the on-shell 
spectral functions in the vacuum. As demonstrated in 
Ref. [15], the off-shell dynamics is important for reso- 
nances with a rather long lifetime in the vacuum but 
strongly decreasing lifetime in the nuclear medium (es- 
pecially uj and <j) mesons) and also proves vital for the 
correct description of dilepton decays of p mesons with 
masses close to the two pion decay threshold. For a de- 
tailed description of the off-shell dynamics and the imple- 
mentation of vector-meson modifications in the medium, 
we refer the reader to Rcfs. [7, 14, 15, 20, 28]. 

In Ref. [7], the PHSD has been extended to include 
hadronic sources for dilepton production from secondary 
multi-meson interactions through the channels -klj — > 
l + l~, ttoi — > l + l~, and pp — > l + l~ . These so-called 
'47T channels' for dilepton production are incorporated in 
the PHSD on a microscopic level rather than assuming 
thermal dilepton production rates and a parametrization 
for the inverse reaction p + + pT — > At: 1 s is incorporated 
by employing the detailed balance as in Refs. [29, 30]. 
By studying the electromagnetic emissivity (in the dilep- 
ton channel) of the hot hadron gas, it was shown in 
Refs. [2, 31] that the dominating hadronic reactions con- 
tributing to the dilepton yield at the invariant masses 
above the <f> peak are the two-body reactions of n + p, 
7r + uj, p + p, and 7T + a\. This conclusion was sup- 
ported by the subsequent study in a hadronic rclativistic 
transport model [32]. Therefore, we implemented the 
above listed two-meson dilepton production channels in 
the PHSD approach in Ref. [7] . In addition, some higher 
vector mesons (p' etc.) were tacitly included by using 
phenomcnological form factors that are adjusted to the 
experimental data. Specifically, we determined the cross 
sections for the mesonic interactions with dileptons in the 
final state using an effective Lagrangian approach follow- 
ing the works of Refs. [2, 32]. In order to fix the form 
factors in the cross sections for dilepton production by 
the interactions of tt + p, it + uj, p + p and rrai, we used 
the measurements in the detailed-balance related chan- 
nels: e + e~ — > 7r + p, e + e~ — > it + ui, e + e~ — > p + p, and 
e + e~ — > TT+a-y. Note that we fitted the form factors while 
taking into account the widths of the p and a\ mesons 
in the final state by convoluting the cross sections with 
the (vacuum) spectral functions of these mesons (using 
the parametrizations of the spectral functions as imple- 
mented in the PHSD). In Fig. 5 of Ref. [7] we presented 
the resulting cross sections that were implemented in the 



PHSD. Contributions of these channels to the dilepton 
invariant mass spectrum in nucleus-nucleus collisions at 
top SPS and RHIC energies and comparisons to available 
data have been reported in Refs. [6, 7]. 



III. THE EXTENDED STATISTICAL 
HADRONIZATION MODEL (SHM) 

Dilepton production has also been studied in the 
framework of the Statistical Hadronization Model (SHM) 
at SPS [33] and RHIC energies [10]. Even in the absence 
of a dynamical evolution of the fireball, both analyses 
find fair agreements with measurements. In the present 
work, we extend the previous studies to LHC energies 
in order to model in detail the contribution of dileptons 
stemming from the so-called "freeze-out cocktail" . These 
results can then be compared with our transport sim- 
ulations in order to quantify the effects from partonic 
degrees of freedom and the dynamical evolution of the 
produced hot and dense matter. 

The SHM is a suitable reference model for obtain- 
ing additional dynamical information about the colli- 
sions besides the particle spectra at chemical freeze-out. 
A corresponding comparison for dilepton spectra from 
Au + Au collisions at the top RHIC energy was presented 
in Ref. [7]. Here we aim at a similar comparison for 
Pb+Pb collisions at LHC energies which requires to spec- 
ify the intrinsic parameters of the statistical hadroniza- 
tion model. Since the central rapidity region is, to a good 
approximation, net charge free in heavy ion collisions at 
the LHC energies, we set all chemical potentials to zero 
in this analysis. Based on findings at SPS and RHIC, we 
expect the 7s (strangeness suppression) factor to be close 
to unity in central Pb + Pb collisions at LHC energies. 
This was recently confirmed in Ref. [34] where measure- 
ments of ratios of hadronic rapidity densities were com- 
pared with statistical model predictions [35] evaluated 
with the thermal parameter T = 164McV and 7s=l. All 
ratios of kaons and multi-strange hyperons to pions were 
correctly predicted, indicating that strangeness is indeed 
chemically equilibrated (i.e. 7s = 1) in central Pb + Pb 
collisions at the LHC energies. Therefore, the thermal 
state of the produced fireball can be specified with two 
free parameters: the temperature T and volume V. We 
do not attempt to extract the chemical freeze-out tem- 
perature from the data but use the predicted asymptotic 
value of T=170 MeV [36] instead. The overall normal- 
ization is determined from the measured transverse mo- 
mentum spectra of various hadrons (see below). 

The SHM - describing a static single fireball - is rather 
reliable in predicting the phase-space integrated relative 
yields of different hadrons. However, since experiments 
measure electrons and positrons at mid-rapidity only and 
above some minimum transverse momentum, we need to 
emulate the initial state dynamics such that the detector 
acceptance is correctly accounted for. 

The ALICE collaboration has measured the transverse 
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momentum spectrum of 7r , K and p in the 0-5% most 
central Pb+Pb collisions at 2.76 TeV [34] while H~ and 
f2~ spectra are measured in the 0-20% most central col- 
lisions. We can estimate the multi-strange hyperon ra- 
pidity densities in the 0-5% most central collisions with 
help of the combined ALICE and CMS measurement of 
charged hadron rapidity density as a function of central- 
ity (see Figure 5 of [37]). 

Fitting the above-mentioned charged hadron data by a 
second-order polynomial, we obtain the following scaling 
factor: 



N, 



5^20% 



dN°r 5% ,dN°- 20% 

cn I cn 



1612/1203= 1.34, 



(1) 

by which we multiply the hyperon transverse momen- 
tum spectra in the 0-20% most central collisions in order 
to estimate the corresponding spectra in the 0-5% most 
central collisions. 

We then assume that the hadron yield in Pb + Pb col- 
lisions arises as a superposition of Np art /2 independent 
nucleon-nuclcon collisions and determine the average el- 
ementary volume V pp = 2.70 fm 3 such that our pion 
yield, which scale like Np aTt V P p/2, agree with the mea- 
sured pion yield in central Pb + Pb collisions. We further 
take Np art =383 independent nucleon-nucleon collisions 
for each 0-5% most central Pb + Pb collision. 

The total mass M c i ust of the hadron cluster in the fire- 
ball can be calculated by adding the energies of individ- 
ual hadrons. From the equipartition theorem for a sys- 
tem in thermal equilibrium that energy is shared equally 
among its degrees of freedom, we can assume that the 
"thermal" momentum of the cluster |j? c iust | equals on av- 
erage to the mass of the cluster. Since the longitudinal 
direction is dominantly governed by the initial state par- 
ton dynamics, we only apply this to the cluster's trans- 
verse momentum, i.e., its average transverse momentum 
is (p^ ust ) = -^-Mciust. In each event, we thus first cal- 
culate Mdust and sample a Gaussian distribution for the 
cluster's transverse momentum such that the mean value 
of the Gaussian equals to ^/2/3M c i ust while the width of 
the distribution is taken as M c i ust / 

For the distribution of the longitudinal momentum of 
the cluster, we rely here on the simple Landau scaling, 
which is in fair agreement with the measurements at 
least up to the RHIC energies, for a discussion, see e.g. 
Ref. [38] . According to the Landau scaling, the width of 
the rapidity distribution of pions can be estimated from 
the simple formula 



<U/i)=ln( 



2m„ 



Trip = 0.94GeV, 



(2) 



leading to <jy (2. 76TeV)=7. 26 which we take as the width 
of the clusters' rapidity distribution in our calculations. 
We note that in our previous work at the RHIC ener- 
gies, we have fitted the width of the clusters' rapidity 
distribution with the value <7^(200GeV)=4.2 which is in 
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FIG. 2: (Color on-line) Transverse momentum spectra of 
ir~ , K~ , 3 _ and Q~ (squares) in Pb + Pb collisions at 2.76 
TeV [34] measured by the ALICE collaboration in comparison 
to our statistical model calculations (lines). 



fair agreement with the value a* (200GeV)=4.7 obtained 
from the Landau scaling. All hadrons are then boosted 
to the laboratory frame where the cluster moves with the 

clust 



momentum p c i ust = (p$ , y c iust \J Pj • + M^ lust ). 

A comparison of our calculated transverse momentum 
spectra of ir~ , K~ , S~ and il~ in central Pb + Pb col- 
lisions with ALICE measurements is presented in Fig- 
ure. 2. We have omitted here the transverse momentum 
spectrum of p as it is not in agreement with that mea- 
sured in experiments. Nevertheless, the midrapidity data 
for the different hadron species is rather well described 
such that we can continue with the study of hadronic 
dilcpton decays. 

The form factors for dilcpton decays of the light 
hadrons employed here are the same as in our previ- 
ous study [10]. Since the low-mass dilcpton spectrum 
is strongly dominated by hadronic decay channels, the 
SHM already allows us to evaluate the e + e~ pair spec- 
trum from the different hadronic decays, i.e. the usual 
'hadronic cocktail'. On the other hand, above invariant 
masses of about 1 GeV the dominant hadronic channels 
for dileptons are the correlated and un-corrclated semi- 
leptonic decays of open charm and open beauty mesons, 
whose relative populations are modeled within the ex- 
tended SHM while the total open charm and beauty cross 
sections are taken from elementary p+p collisions at the 
same bombarding energy and scaled with the number of 
binary nucleon-nucleon collisions (see below). 

The uncertainty in the total charm and beauty produc- 
tion cross sections, especially at LHC energies, is sizable. 
In order to estimate these cross sections, we use differ- 
ent (but related) models to achieve a good reproduction 
of the data from pp collisions at v / s NN =2.76 TeV as de- 
scribed in the following Section IV. The resulting cross 
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sections are used to calculate the production of D- and 
D -mesons in Pb-\- Pb collisions at -\/^nn — 2.76 TcV. Con- 
tributions from lcptonic decays of these heavy mesons are 
then added to the dilepton spectra obtained from SHM 
and PHSD described in Section VI. 



A. 



IV. OPEN CHARM AND BEAUTY 



Open charm and beauty in proton-proton 
collisions 



There are experimental indications [39, 40] that heavy 
flavor production might also exhibit statistical features 
in high energy nuclear as well as in e + e~ collisions, at 
least for the relative abundances of open charm [41] and 
possibly also for that of open beauty mesons [42]. Since 
heavy-flavor production occurs early in the reactions and 
is not likely formed thermally, we thus extend the sta- 
tistical hadronization model to include open charm and 
beauty mesons such that their relative yields are eval- 
uated with the SHM (for some freeze-out temperature 
T) while their total numbers are estimated from mea- 
sured open and beauty production cross sections. For 
the D s and B s states, their production cross sections 
are multiplied by the factor 7^ ard =0.3 relative to those 
for nonstrange heavy mesons to take into account the 
strangeness suppression in hard proton-proton scatter- 
ing. Furthermore, we take into account only the 6 lowest 
mass states of open charm mesons in the extended SHM 
since they already account for more than 90% of the total 
open charm production cross section and also not much 
are known about the higher mass excited states. 

With the above extended SHM, we study the open 
charm and beauty production for the two freeze-out tem- 
peratures T=170 MeV and 150 McV while the chemical 
potentials are kept zero. The resulting relative normal- 
ized primary and final (after strong decays) production 
probabilities for the open charm states (calculated within 
the SHM) are given in Table I. The ratios of final open 
charm hadrons can then be calculated and compared with 
those measured by the LHCb and ALICE Collaborations 
in p + p collisions at •\Ann = ? TeV [43-45] as shown in 
Table II. As one can see, all ratios are reproduced (within 
errors) for both temperatures of 150 and 170 MeV. We 
choose here, as for the light hadrons, T=170 McV as 
our temperature for all hadron species considered in the 
analysis. We note that the feed-down contribution from 
beauty mesons was neglected in the above analysis since 
this contribution is estimated to be of the order of 10% 
of the total charm production cross section [46] . 

In the above, we have taken the total open charm pro- 
duction cross section in p + p reactions at 7 TcV to be 



7ToV 



6.4 mb 



from fitting a Gaussian to the combined measurements 
of ALICE and LHCb D on meson rapidity distributions. 
Again, the longitudinal dynamics of open heavy flavor 



Hadron 


primary 


final state 




T=170 MeV 


T=150 MeV 


T=170 MeV 


T=150 MeV 


D+ 


0.186 


0.200 


0.273 


0.286 


D° 


0.190 


0.204 


0.646 


0.642 


D s 


0.0333 


0.0331 


0.0810 


0.0751 


D* + 


0.269 


0.258 








D"° 


0.274 


0.263 








D* 3 


0.0474 


0.0420 









TABLE I: Relative production probabilities of open charm 
mesons before (primary) and after (final) strong decays. 



D 



I 

D. 
D* + 



D_ 

D, 



LHCb [43, 44] 



2.20 ± 0.48 
2.07 ± 0.37 
7.67 ± 1.67 
0.94 ± 0.22 
3.48 ± 0.93 
3.70 ± 0.84 



ALICE [45] 



2.09 
2.08 

1.00 



T=170 MeV 



2.40 
2.37 
7.98 
0.99 
3.32 
3.37 



T=150 MeV 



2.49 
2.25 
8.55 
0.90 
3.44 
3.81 



TABLE II: Ratios of measured [43-45] open charm yields in 
p + p collisions at 7 TeV compared with SHM model calcula- 
tions for the temperatures T=170 MeV and T=150 MeV. 



mesons are evaluated in the model of Rcf. [10] while the 
transverse momentum distributions are predicted by the 
MC@sHQ model [47-49], explained in Section IV B. Our 
estimate for the total charm production cross section (6.4 
mb) is comparable to the preliminary LHCb estimate of 
6.1 ± 0.93 mb obtained by tuning PYTHIA to the LHCb 
measurement and then integrating over the whole phase- 
space [44] . An extrapolation of the ALICE collaboration 
measurement at mid-rapidity to full phase-space yields, 
however, a cross section of CTgJ eV = 8.5+ 2 ; 4 mb [50]. The 
ATLAS collaboration gives, on the other hand, for the to- 
tal cross section a preliminary value of 7. 13^2 ° m t> [51]- 
The three LHC experiments agree within the uncertain- 
ties. Pcrturbativc QCD calculations were also found to 
agree with the data and with our estimate for the total 
charm production cross section in p+p reactions at 7 TeV 
(see [45, 52] and references therein). 

For the total charm production cross section at 2.76 
TeV reads, it is estimated to be ct? c 76TcV =4.8±0.8 mb 
by the ALICE Collaboration. Using the ratio of the two 
ALICE cross sections (at 7 TeV and at 2.76 TeV), we 
find that our estimate for the total cross section of charm 
production hi p + p collisions at the lower LHC energy is 



a c 2 c 76TcV = 3.6 mb, 

which is used in our study for heavy ion collisions at 
2.76 TeV. 

For the open charm decays, we have implemented in 
detail the form factors for each of the semi-leptonic decay 
channels that are inlcuded in our analysis. A discussion 
and presentation of relevant formulae for the open charm 
form factors can be found in the Appendix A. For beauty 
decays, no form factors were employed. 
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The properties of open beauty mesons are not well 
known experimentally and thus we use the known prop- 
erties from the PDG [53] and take the missing informa- 
tion from model calculations. In the case of open beauty, 
we take into account 18 lightest states and their anti- 
particles. Only the lowest lying states B + , _B°, B s and 
B* and their anti-particlcs decay semi-lcptonically into 
electrons, while the excited states decay under the strong 
interaction into the lower mass states. It is important to 
take into account also the excited states in the analysis in 
order to properly evaluate the total beauty cross section 
and the relative abundances of the lowest lying B meson 
states. We ignore in this work the 6-baryons whose pro- 
duction cross section is expected to be small (as 10% of 
°bb) [53] and thus our a^ h refers to the total cross section 
for the production of open beauty mesons. 

The properties (masses and decay channels) as well 
as the references to the data and models that we use 
for our open beauty calculations are listed in Table VI 
in Appendix A. In all cases we have taken a weighted 
average over the different estimates as input to our study, 
and a recent compilation for the numerical values of the 
open beauty masses can be found in [54] . 

Measurements of the total beauty cross section is not 
yet available at LHC energies, hence we relate it to the 
charm cross section and use available theoretical esti- 
mates for the ratio of these cross sections to estimate a^ h 
in p+p collisions at 2.76 TeV. We note that perturbative 
QCD predictions for this ratio suffer from large uncer- 
tainties [55]. Using here the same model and parameters 
as those applied for the calculation of the energy loss of 
heavy quarks in the next Section, we obtain the ratio 
cr E c/ob b =40 in central Pb + Pb collisions at 2.76 TeV. 



B. Energy loss of D- and B-mesons in Pb + Pb 
collisions 

The total cross sections and spectra of D-mesons ob- 
tained in the previous Section from the available exper- 
imental data are for elementary p+p collisions. We are 
interested in the contribution of the decays of correlated 
and uncorrelated D + D and B + B pairs to the dilepton 
spectra in heavy ion collisions, but simple binary scaling 
of the spectra from p+p collisions is expected to overesti- 
mate their contributions. This is because of the neglect of 
the loss of correlation and the modification of the spectra 
of heavy hadron pairs due to the partonic energy loss [49] 
and hadronic rescattering [56]. To take into account the 
collisional and radiative energy loss of charm and beauty 
quarks in central Pb+Pb collisions, we use the event gen- 
erator for heavy quarks and mesons of P.B. Gossiaux et 
al. - Monte Carlo for Heavy Quarks (MC@sHQ). This 
model - described in detail in Refs. [47-49] - assumes 
that in heavy-ion collisions the initial momentum distri- 
bution of charm and bottom quarks is identical to that 
in p+p collisions which arc predicted by the FONLL 
PQCD calculations [57] and verified experimentally by 



the STAR and PHENIX collaborations [58]. With the 
expansion of the produced quark-gluon plasma described 
by the hydrodynamic model, the interaction of heavy 
quarks with the plasma constituents is then calculated 
from the PQCD one-gluon-exchange diagrams but with 
a running coupling constant and an infrared regulator ad- 
justed to reproduce the same energy loss as in the hard 
thermal loop calculations. Both collisional and radiative 
energy loss are taken into account and with the latter cor- 
rected by the Landau-Pomeranschuk-Migdal effect. This 
model has been shown to describe very well the nuclear 
modification factor Raa and the elliptic flow vi of heavy 
mesons (non-photonic electrons) measured at RHIC and 
more recently at LHC by the ALICE Collaboration [59] . 

Here, we employ the transverse momentum distribu- 
tions of the D- and B-mesons produced in central Pb+Pb 
collisions as calculated in the MC@sHQ event generator 
as an input for our calculation of the dilepton yield by 
the semi-leptonic decays of D + D and B + B meson pairs. 
In particular, the decay of a single D (D) meson is sim- 
ulated and the produced electron (positron) is combined 
with the positron (electron) produced in the decay of the 
D (D) meson coming from the same production vertex to 
form dilepton pairs from the correlated D + D decays. A 
-D-meson-decay electron is combined with a positron pro- 
duced in the decay of a D-meson coming from a (random) 
different production vertex but from the same event to 
form dilepton pairs from the uncorrelated D + D decays. 
Furthermore, our prescription for estimating the loss of 
correlations due to the interactions of charm in the fire- 
ball is to multiply the yield from the correlated pairs by 
the suppression factor R?aai where Raa is the measured 
nuclear modification factor for the D-mesons. A similar 
procedure is used for treating the semi-lepton decays of 
£?-mcsons. 



V. CHARMONIA 
A. Charmonia in proton-proton collisions 

Quarkonia are important (and experimentally clean) 
sources of dilepton in the high mass region and can also 
be used to constrain and cross-check our results regarding 
open charm production. Unfortunately, the SHM is not a 
suitable model for studying the production of the hidden 
charm hadrons. Therefore, we include the 3/ip and "0' 
mesons in our analysis as an external input taken from 
the experiments, i.e. we do minimal modeling for these 
hadrons and take basically all necessary information from 
available data. 

The ALICE collaboration has measured the 3/tp ra- 
pidity and transverse momentum distributions in p + p 
collisions at 2.76 TeV [60]. We fit a Gaussian to the mea- 
sured rapidity distribution in order to extract the inclu- 
sive J/ip production cross section of 29 fib for this collision 
system as shown by the solid line in Fig. 3. According 
to the color evaporation model, the relative production 
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p+p, s =2.76 TeV, lyl<0.5 



a= 29.32 lib 




-8 -6 -4 -2 2 4 6 8 



FIG. 3: (Color on-line) 3 /if) rapidity distribution in p+p colli- 
sions at 2.76 TeV measured by the ALICE Collaboration [60] 
(full spheres) and reflection of the data (open spheres). Solid 
line is a Gaussian fit to the data based on a total J /if) pro- 
duction cross section of 29.3/^6 for this collision system. Also 
shown by squares is our simulated 3 /if) rapidity distribution 
as described in the text. 

probabilities of different quarkonia are beam-energy inde- 
pendent at ultra-relativistic energies (for a recent review, 
see Ref. [61]) and, once the 3 /if) cross section is known, 
we can estimate the production cross section of if)' in 
p + p collisions at 2.76 TeV with the help of the world 
average [62] 2^^=8.1%. 

The transverse momentum spectrum of 3/ipcan also 
be determined from the ALICE data [60] , and we use the 
same spectrum for J /if) and if)' both in p + p as well as 
in Pb + Pb collisions. The parameters characterizing the 
transverse momentum spectra of hidden and open heavy 
flavor hadrons are collected in Table V in Appendix B. 
We mention here that the details of the pr spectrum 
play only a minor role in di-electron radiation from the 
quarkonia, while in the case of D and B meson decays 
the results are much more sensitive to the transverse mo- 
mentum spectrum. Due to this reason, we can use the 
same px profile for quarkonia in p + p and A + A while 
for open heavy flavor decays the two systems have to be 
modeled separately. 

The longitudinal momentum distributions of all heavy 
flavor mesons are modeled phcnomenologically as ex- 
plained in detail in our earlier work [10]. In the present 
calculations, we use the same relevant parameters (the 
choice of parton distributions, which determines the 
width of the clusters' distribution in longitudinal direc- 
tion) which we have fixed at RHIC energies. Further- 
more, we use the CT10 [63] parametrization for the par- 
ton distributions in all calculations. The calculated ra- 
pidity distribution of 3 ftp's within this model (with the 
prc-dctcrmincd cross section) is compared to the ALICE 
data and to the Gaussian paramcterizations in Fig. 3. We 



find that the measured J /if) rapidity distribution is nicely 
reproduced and will use the same procedure to evaluate 
the longitudinal momentum distributions of all heavy fla- 
vor hadrons in this work. We mention here (without 
explicit presentation) that the measured rapidity distri- 
butions of open charm hadrons at both mid-rapidity [45] 
and forward rapidities [43, 44] in p + p collisions at 7 
TeV are also in good agreement with the results from 
this model. 



B. Nuclear modification of charmonia 

To estimate the J/4" contribution to the dilepton spec- 
trum in heavy ion collisions, we need information on the 
J /if) nuclear modification factor in these collisions. This 
can be determined using the approach of Refs. [64, 65] 
that includes charmonium production from both initial 
hard nucleon-nucleon scattering and regeneration from 
charm and anticharm quarks in the produced QGP. For 
the initially produced J/^f, their number is proportional 
to the number of binary collisions between nucleons in the 
two colliding nuclei. Whether these J/^'s can survive af- 
ter the collision depends on many effects from both the 
initial cold nuclear matter and the final hot partonic and 
hadronic matters. The cold nuclear matter effects include 
the Cronin effect of gluon-nucleon scattering before the 
production of the primordial J /"J from the gluon-gluon 
fusion [66], the shadowing effect due to the modification 
of the gluon distribution in a heavy nucleus [67] and the 
nuclear absorption by the passing nucleons [68-70]. For 
the hot partonic and hadronic matter effects, they in- 
clude the dissociation of charmonia in the QGP at tem- 
peratures higher than the dissociation temperature and 
the thermal decay of surviving charmonia through inter- 
actions with thermal partons in the expanding QGP. In 
studying these effects, medium effects on the properties 
of the charmonia and their dissociation cross sections are 
further taken into account by using the screened Cornell 
potential model [71] and the NLO pQCD [72]. This leads 
to a dissociation temperature of ~ 300 MeV for J /if) and 
~ T c = 175 MeV for where T c is the QCD phase 
transition temperature. Although the J/\& can survive 
above Tc, its thermal decay width is not small with a 
value of - 10 MeV at T = 200 MeV and increasing 
to - 100 MeV at T = 280 MeV. Because of the ap- 
preciable number of charm quarks produced in relativis- 
tic heavy-ion collisions, charmonia can also be regener- 
ated from charm and anticharm quarks in the QGP. A 
rate equation is then used to include the effect of ther- 
mal dissociation and regeneration of charmonia. Since 
charm quarks are not expected to be completely ther- 
malized - neither chemically nor kinetically - during the 
expansion of the hot dense matter, a fugacity parameter 
and a relaxation factor are introduced to describe their 
distributions. Modeling the evolution of the hot dense 
matter produced in relativistic heavy-ion collisions by a 
schematic viscous hydrodynamics [64], the average tern- 
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FIG. 4: (Color on-line) Dielectron invariant mass spectrum 
from various different sources evaluated for p + p collisions at 
2.76 TeV within the SHM for the invariant mass range up to 
1.2 GeV (upper window) and 4 GeV (lower window). 
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FIG. 5: (Color on-line) Same as Fig.4 from the PHSD ap- 
proach. 



perature of the initially produced QGP in central Pb+Pb 
collisions at ^/snn = 2.76 TeV is ~ 311 MeV if the ini- 
tial thermalization time is taken to be 1.05 fm/c. The 
resulting nuclear modification factors for J/^>(1S) and 
\I/'(2S) are 0.30 and 0.39, respectively, if the shadowing 
effect is included. These values increase to 0.57 and 0.95, 
respectively, in the absence of the shadowing effect. 



VI. RESULTS 

A. The baseline: p+p collisions 

As a baseline for our dilepton studies, we first show in 
Fig. 4 the SHM results for the invariant mass spectrum 
of e + e~ pairs produced in p + p collisions at -\/snn = 



2.76 TeV. Qualitatively, the spectra are very similar to 
those from the top RHIC energy of V«nn= 200 GcV I 10 ] 
but enhanced by roughly a factor of 2. These results 
have been obtained with a transverse momentum cut of 
150 MeV and a rapidity cut of \y e \ < 0.88 for the sin- 
gle electrons, similar to the nominal acceptance cuts of 
the ALICE detector [73]. Also, the centrality and pr de- 
pendent paramctcrizations of the ALICE detector resolu- 
tion [73] have been used for determining the mass resolu- 
tion. We have further included the contribution from the 
initial state Drell-Yan process in the collinear factoriza- 
tion approach of the next-to-leading-order perturbative 
QCD [74-76] (see also Ref. [77] for a recent discussion 
at LHC energies), which is expected to be compatible 
with the experimental data within a factor of two. We 
find that the Drell-Yan process contributes less than 1% 
in the dilepton invariant mass range considered in our 
studies and thus discard its contribution. 
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TABLE III: Meson decay channels for di-electron production. 



Next, we present in Fig. 5 the PHSD results for the 
invariant mass spectrum of e + e~ pairs produced in p + p 
collisions at v / *nn = 2-76 TeV, employing the same trans- 
verse momentum and rapidity cuts for the single elec- 
trons, i.e. pfp > 150 MeV and \y e \ < 0.88. We observe 
generally a good agreement between the two models. The 
minor differences in the line shapes stem from i) the 
production of hadrons in proton-proton collisions in the 
PHSD from string decays versus the thermal distribution 
of primary produced particles (combined with the feed- 
down from resonances) in the SHM and ii) the slightly 
different compositions of dilepton-decay channels taken 
into account in the SHM and the PHSD. The list of the 
channels for the meson decays to dileptons - contributing 
to the di-electron spectrum in proton-proton collisions - 
is presented in Table III. The (non-dominant) channels, 
which are taken into account in the SHM but not in the 
PHSD, are marked with asterisks. 



B. Results for central Pb + Pb collisions at 
v/I NN =2.76 TeV 

We first present the results for dilepton production 
from the SHM at y / s NN =2.76 TeV and then compare to 
the spectra from the PHSD transport approach, which 
incorporates the same cross sections for cc as well as bb 
pair production as specified above, but also allows for 
an estimate of hadronic interaction channels as well as 
in-medium spectral functions and the radiation from the 
sQGP. We concentrate on most central Pb + Pb colli- 
sions, since the uncertainties in the number of partici- 
pants Np art are low and hadronic in-medium effects are 
most pronounced. 

In Fig. 6 we display our results from the SHM for the 
dilepton low mass sector (upper window) as well as for 
the invariant mass range upto 4 GeV (lower window). 
We recall that the SHM does not include any dynamical 
effects and only accounts for the final hadronic 'cocktail'. 
The results for di-electron decays of light mesons are ob- 
tained by a scaling of p+p results with the number of par- 
ticipants as measured by the ALICE Collaboration [23]. 
The only difference between the statistical hadronization 
model calculations for the elementary and the heavy-ion 
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FIG. 6: (Color on-line) Dielectron invariant mass spectrum 
from various sources evaluated for the 0-5% most central Pb+ 
Pb collisions at ^/s]7n = 2.76 TeV within the SHM for the 
dielectron invariant mass range up to 1.2 GeV (upper window) 
and 4 GeV (lower window). Only the un-correlated D and B 
meson contributions are shown since almost all correlations 
are destroyed in central Pb + Pb collisions. 



reactions is that 7s=0.6 was used for p+p collisions, while 
7s=0.95 was used in Pb + Pb collisions, because the lat- 
ter value better fits the transverse momentum spectrum 
of light mesons and multistrange hyperons produced in 
Pb + Pb collisions at Vs NN =2.76 TeV. 

We point out explicitly that rescattcring is fully in- 
corporated in our calculations for the contributions from 
open charm and beauty decays as well as in the dilep- 
ton decays of charmonia (in Fig. 6 and in all subsequent 
figures in this paper). While heavy- flavor mesons are pro- 
duced fully back-to-back in p + p collisions (cf. Fig. 4), 
we assume that these initial correlations are washed out 
by final-state interactions in central Pb + Pb collisions 
since the survival probability of the D and D correlation 
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FIG. 7: (Color on-line) Same as Fig.6 from the PSHD with 
vacuum spectral functions for the vector mesons. 



is approximately given by Raa{D) 2 , where Raa is the 
-D-meson nuclear modification factor. Also, the energy 
loss of the produced D and B mesons was taken into 
account, both in the correlated and uncorrelated contri- 
butions, according to the Nantes model as described in 
Section IV B. Additionally, the nuclear modification of 
charmonia was taken into account by using the predic- 
tions from Section VB. 

Results from the static calculations within the SHM 
can be confronted with those from the dynamical trans- 
port approach PHSD that are displayed in Fig. 7. In 
general, the low mass spectrum from the PHSD shows 
a similar channel decomposition as in the SHM, while 
the intermediate mass region is enhanced substantially 
by the partonic radiation channels which are dominated 
by the quark-antiquark annihilation channel. We note 
that although the contribution of the dileptons from the 
correlated decays of D + D and B + B mesons also peaks 
at the intermediate mass region, its total yield is low. 
This is due to the fact that the (partonic) rescattering 




M [GeV/c 2 ] 

FIG. 8: (Color on-line) Upper part: Dielectron low invariant 
mass spectrum from various different sources evaluated for 
10% central Pb+Pb collisions at % /s NN =2.76 TeV within the 
PHSD approach including in-medium spectral functions for 
the vector mesons. The lower part displays our results in the 
mass range up to 4 GeV. 



of the charm and bottom quarks in central Pb + Pb col- 
lisions practically destroys all the correlation of the ini- 
tially produced heavy-quark pairs in p + p collisions. It 
is encouraging to see that (Fig. 7) the sQGP radiation is 
clearly visible in the intermediate mass region, and this 
provides the possibility to measure the dilepton radia- 
tion from the sQGP, making the study of its properties 
experimentally accessible. A similar conclusion was also 
made in Ref. [7] for heavy ion collisions at the top RHIC 
energy. 

Another interesting question is the modification of the 
vector-meson spectral functions by hadronic in-medium 
effects which was found to be essential for describing the 
NA60 dilepton data in the low mass regime from heavy 
ion collisions at 158 A-GeV available from the SPS [6]. 
In Fig. 8 we present the dilepton invariant mass spec- 
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FIG. 9: (Color on-line) Dielectron low invariant mass spec- 
trum from various different sources evaluated for most cen- 
tral Pb+Pb collisions at ^/s NN =2.76 TeV within the PHSD 
approach with vacuum as well as with in-medium spectral 
functions for the vector mesons. 



trum from the PHSD in the collisional broadening sce- 
nario for the most central Pb+Pb collisions at y / i NN = 
2.76 TeV. It is seen that the spectra in the intermediate 
mass region are not much affected by the medium mod- 
ified vector meson spectral functions due to the domi- 
nance of partonic contributions and correlated dilcptons 
from Z?-meson decays. A direct comparison of the two 
above calculations is displayed in Fig. 9, where we focus 
on the low mass sector. A moderate in-medium effect 
(within 20% of the total dilepton yield) in the mass win- 
dow from 0.3 to 0.6 GeV due to the modification of the 
vector mesons can be identified. 

We next study the effect of a transverse momentum 
cut on the channel decomposition of the dilepton spec- 
trum. Shown in Fig. 10 are the PHSD results for the 
dilepton yield dN/ dM in very central Pb + Pb collisions 
at •y / s NN =2.76 TeV and mid-rapidity (as in Fig. 7) but 
with the additional cut on the transverse momentum of 
single electrons of pt > 1 GeV. We find that the sQGP 
radiation signal in this case is even more pronounced due 
to the suppression of the heavy meson contribution to 
the dilepton yield in the intermediate masses AI = 1 — 3 
GeV. Thus the pt cut is beneficial to the extraction of 
the sQGP signal in the intermediate mass region. 

On the other hand, we see in the upper part of Fig. 10 
that the low-mass part of the spectrum is strongly mod- 
ified by the cut: the uj and <f> peaks become more promi- 
nent while the total yield is strongly suppressed. This 
loss of statistics might impair the detailed measurements 
of the vector meson spectral shape and of the relative 
meson yields. Thus we suggest that, although the cut 
on px is profitable for the measurement of the intermedi- 
ate mass dilcptons in order to study the sQGP radiation, 
the low-mass spectrum of the dileptons (M < 1.1 GeV) 
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FIG. 10: (Color on-line) Upper part: Dielectron low invariant 
mass spectrum from various different sources evaluated for 
the central Pb+Pb collisions at ^/s NN =2.76 TeV within the 
PHSD approach including vacuum spectral functions for the 
vector mesons with the transverse momentum cut pr > 1 GeV 
on the single electron. The lower part displays our results in 
the mass range up to 4 GeV. 

should better be measured down to lowest possible values 
of px- 



VII. SUMMARY 

In this study, we have extended our previous investi- 
gations of dilepton production at top SPS and top RHIC 
energies [6, 7, 10] to p + p as well as central Pb + Pb 
collisions at A /s NN =2.76 TeV available at the LHC by 
employing an extended statistical hadronization model 
(SHM) and the dynamical PHSD transport approach. 
The SHM serves to independently determine the back- 
ground from the 'hadronic cocktail', while the PHSD ad- 
ditionally provides information on dynamical contribu- 
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tions to the dilepton spectrum, i.e. a low-mass enhance- 
ment due to vector meson modification by hadronic scat- 
terings as well the contribution from partonic reactions. 
Both of the above dynamical effects are not visible in the 
final hadronic decay contributions. 

In extending our previous investigations, a few new 
developments have been incorporated apart from a read- 
justment of the two SHM parameters (T and V) to LHC 
energies: i) an inclusion of bottom quarks and antiquarks 
as well as the corresponding bottom mesons and their 
leptonic decay modes, ii) an appropriate model for the 
energy loss of open charm mesons in central Pb + Pb col- 
lisions [47-49] , and iii) a dynamical approach for the char- 
monium suppression and recreation in the partonic and 
hadronic phases [64, 65] that supersedes the earlier stud- 
ies within the HSD approach at lower energies [78] . With 
respect to the partonic channels for dilepton production 
in the PHSD no extensions had to be incorporated as 
well as with respect to the overall reaction dynamics of 
Pb + Pb collisions in the transport approach. 

We find a reasonable agreement between the differ- 
ential mass spectra of dileptons from hadronic decays 
calculated in the PHSD approach and the SHM model. 
Dynamical contributions or medium effects due to a 
broadened p spectral function is seen to give a rela- 
tively small correction in the low-mass dilepton sector 
from 0.3 — 0.6 GeV, since the hadronic final state in- 
teractions are less important at LHC energies. How- 
ever, pronounced traces of the partonic degrees of free- 
dom are found in the PHSD results for the intermedi- 



ate mass regime, i.e. between the (j> and J/\& decay 
peaks, with the dilepton yield from the strongly inter- 
acting quark-gluon plasma (sQGP) even exceeding that 
from the semi-leptonic decays of open charm and bottom 
mesons. Additionally, we observe that a transverse mo- 
mentum cut of 1 GeV/c for the leptons further suppresses 
the relative contribution of the heavy meson decays to the 
dilepton yield in the intermediate mass region, such that 
the sQGP radiation strongly dominates the spectrum for 
masses from 1 to 3 GeV, allowing for a closer look at the 
electromagnetic emissivity of the partonic plasma in the 
early phase of Pb+Pb collisions. 
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Appendix A: D meson form factors 

The form factors for D° — > K~e + v e decays have been 
measured by CLEO [79], FOCUS [80], BELLE [81] and 
BaBar [82] Collaborations. These can be conveniently 
parametrized in the modified pole ansatz [83] 



U(Q 2 ) = 



/+(0) 



(l - 



(Al) 



We take a weighted average a=0.35 of the measured [79- 
82] a and tod* =2.112 GeV as our input for the calcula- 
tions and use the same parameters for the corresponding 
D ± decays. 

The measured [79-81] form factors for D —> irlv e de- 
cays can be parametrized in the simple pole approxima- 
tion 



/ + (? 2 ) 



/+(0) 



(1- 



(A2) 



in which we take a weighted average m po i e =1.92 GeV as 
input for our calculations. 

The other D — > Pev e form factors have not been mea- 
sured yet and thus we will rely on theoretical modeling 
here. We use the formalism [84] introduced by Fajfer 
and Kamenik in which the form factor for the decay of 
a D meson (H) has the same form as Eq. (Al) but the 
parameters of the distribution are fixed to the masses of 
the known (nearest) resonances. Explicitly, it is given by 



/+(0) 



(1 



-)(1 



-) 



(A3) 



where a = m 2 H , /m 2 ,,. The nearest pole masses ran* 
and m H u that fix the value of the a parameter are 
process dependent and are listed in the Table I of [84]. 
The predictions [84] of the Fajfer-Kamenik model for the 
scmi-lcptonic branching fractions D + — > {rj' , rf\e + v e were 
recently confirmed by the measurements [85, 86] of the 
CLEO Collaboration. With the help of the form factors, 
the q 2 distribution of the D — > pseudo-scalar + e + v e 
can be calculated from 



^-\px(l 2 )\ 3 \U(<l 2 )\ 2 , 



(A4) 



where X = {K, ir, 77, 77', /o}. The parameters we have 
used for the D — > Pev e transitions are listed in Table IV. 

The -D(s) — > Vev with V denoting a vector meson (A"* 
or <fi) decays can be parametrized with the help of three 
form factors [88, 891 



V(q 2 ) 



V(0) 



(1 — x)(l — ax) 

MO) 



1 



M<f) 



A 2 (0) 



(A5) 



channel 


m po ie [GeV] 


Q 


references 


D — > nev e 


1.92 





[79-81] 


D — > Keu e 


2.112 


0.385 


[79-82] 


D — > r\ ev e 


2.010 


0.74 


[84] 


D — > r/'ev,. 


2.010 


0.74 


[84] 


D a —¥ rjeu e 


2.112 


0.75 


[84] 


D s — > r)'ev e 


2.112 


0.75 


[84] 


D s Kev e 


2.010 


0.74 


[84] 


D a / (980)ei/ e 


1.7 





[87] 



TABLE IV: List of parameters characterizing the D 
transitions. 



Peu e 



in the helicity amplitude formalism. In the above, x = 
1 2 / m poie with mpoie being the (process dependent) pole 
meson mass and the free parameters of the model (V(0), 
Ai(Q), ^.2(0), a and b) were determined in [89] (see Table 
I of [89] for the list of m po i e and Table II for the values of 
the fit parameters), and we use these parameters in our 
calculations. 

The g-distribution of the total decay rate is evaluated 
with 

dr = d(r + + r_ + r ) 

dy 



dy 

y\Pv(y)\(\H + (y)\ 2 + \H-(v)\ 2 + \Ho(y)\ 2 ) 



where y = q 2 /m 2 H and 



2 \2 



\pv(y)\ 2 = 



(m- H (l - y) + m v ) 
4m 2 H 



(A6) 



(A7) 



is the three momentum of the vector meson in the rest 
frame of the D meson, and the helicity amplitudes read 



H±(y) 

H (y) 



(m H + m v )A 1 (m H y) =p 
niH + my 



2 ^^2m H \ P y(y)\ vimly)> 



nriH + mv 

2 \ A /_2 



2m H m Vy /y 

2m H \pv{y)\ 
m v {m H +m v )y/y 



(m H (l - y) - m v )A 1 (m H y) 



A 2 {m%y). 



(A8) 



We have implemented in detail the form factors for 
each and every of the known open charm semi-lcptonic 
decay channels. The only decay channel that is consid- 
ered isotropic in the momentum space is the D — ¥ p°ev e . 
Due to the large width of the p° , above formalism is not 
applicable for this channel. One could implement the 
ISGW2 paramctrization, but this channel contributes 3% 
of the total semi-leptonic decay width for D° and even 
less for the D + . Omission of this form factor should not 
cause much effect. 



Appendix B: Characteristics of the transverse 
momentum spectra 

We do not calculate the transverse momentum spec- 
trum of heavy hadrons within the SHM but fit a power- 
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Hadron 


system 


Ref. 


B [GeV] 


n 


p T [GeV] 


J/V> 


p+p 2.76 TeV 


[60] 


4.53 


4.61 


[0:8] 


D 


p+p 2.76 TeV 


[47-49] 


2.37 


3.09 


[0:10] 


B 


p+p 2.76 TeV 


[47-49] 


6.38 


3.10 


[0:10] 


D 


Pb+Pb 2.76 TeV 


[47-49] 


2.10 


3.78 


[0:10] 


B 


Pb+Pb 2.76 TeV 


[47-49] 


4.19 


3.24 


[0:10] 



TABLE V: Parameters characterizing the heavy flavor trans- 
verse momentum spectra around mid-rapidity in p+p and 
Pb + Pb collisions at 2.76 TeV. 



Hadron 


mass 


relevant branching fractions 




[53] 


e^VeXc 10.2% 


7-)0 

B' 


[53] 


e + u e X c 10.1% 


T)* + 


[53] 


B 7 997o B e e 0.57o [90] 


B*° 


[53] 


B°7 99% B°e + e- 0.5% [90] 


Bt 


[91-94] 


Bit 


B° 


[91-94] 


Bit 


Bt 


[53] 


B*ir [95] 


B° 


[53] 


B*tt [95] 


B*+ 


[91-94] 


B*tt [95] 


Bf 


[91-94] 


B*tt [95] 


b; + 




53] 


Bn 40% B*n 60% [95] 


Bf 




53] 


Bit 40% B*n 60% [95] 


B s 




53] 


charm 


B* a 




53] 


Bh 


B s o 


[91-94, 96] 


BK[Q7] 


B s i 


[53] 


B*K[97] 




[91-94, 96] 


B*K[97] 


£>s2 


[53] 


Btt 60% B*ir 40% 



TABLE VI: References to the works from which we have taken 
the B meson masses and relevant branching ratios. Unless ex- 
plicitly stated otherwise, the branching fractions to the listed 
channels are 100% 



law function to the available experimental data (J/VO 
and theory calculations (D and B mesons) and use these 
paramctcrizations as our input for the simulations. We 
list here all the relevant parameters related to transverse 
momentum spectra of various hadrons. Explicitly, 



d 2 AT 
dprAy 



C(y) 



Pt 



Pt 



(Bl) 



where C(y) is some rapidity dependent normalization fac- 
tor and the parameters B and n are listed in Table V 
along with the references to the input data. 

We use the same px spectrum for -0' as that for J/ip. In 
contrast to the dilcpton emission from D and B mesons, 
our results do not depend strongly on the details of the 
J ftp transverse momentum spectrum and hence we use 
the same parameters for Pb + Pb collisions as in p + p 
collisions. 

Appendix C: Properties of open beauty mesons 



Properties of open beauty mesons are given in Ta- 
ble VI. 



